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Recent experiments performed on a variety of soft glassy materials have demonstrated that any 
imposed shear flow serves to simultaneously fluidize these systems in all spatial directions [Ovarlez 
et al. (2010)]. When probed with a second shear flow, the viscous response of the experimental 
system is determined by the rate of the primary, fluidizing flow. Motivated by these findings, we 
employ a recently developed schematic mode-coupling theory [Brader et al. (2009)] to investigate 
the three dimensional flow of a colloidal glass, subject to a combination of simple shear and uniaxial 
compression. Despite differences in the specific choice of superposed flow, the flow curves obtained 
show good qualitative agreement with the experimental findings and recover the observed power law 



fSJ describing the decay of the scaled viscosity as a function of the dominant rate. We then proceed 

y—< to perform a more formal analysis of our constitutive equation for different kind of 'mixed' ffows 

consisting of a dominant primary ffow subject to a weaker perturbing flow. Our study provides 
04 further evidence that the theory of Brader et al. (2009) reliably describes the dynamic arrest and 

(— ( mechanical fluidization of dense particulate suspensions. 

h— 5 PACS numbers: 47.57.Qk, 82.70.Dd, 83.10.Gr, 83.60.Df 
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have long been a staple of continuum rheology (such as 
the Giesekus or Oldroyd models [Bird et al. (1987), Lar- 
son (1988)]) and enable e.g. normal forces and secondary 
flows to be addressed in realistic curvilinear experimental 
geometries. 

The first steps towards a unified, three dimensional 
description of colloid rheology have been provided by re- 
cent extensions of the quiescent mode-coupling theory 
(MCT) to treat dense dispersions under flow [Brader et 
al. (2008)]. These developments are built upon earlier 
studies focused on simple shear [Brader et al. (2007), 
Fuchs (2009), Fuchs and Gates (2002), Fuchs and Gates 
(2009)] and capture the competition between slow struc- 
tural relaxation and external driving, thus enabling one 
of the most challenging aspects of colloid rheology to 
be addressed: the flow response of dynamically arrested 
glass and gel states. Given the equilibrium static struc- 
ture factor as input (available from cither simulation or 
liquid state theory [Brader (2006)]), the deviatoric stress 
tensor cr{t) may be determined for any given velocity 
gradient tensor K{t). However, implementation of the 
theory has been hindered by the numerical resources re- 
quired to accurately integrate fully anisotropic dynam- 
ics over timescales of physical interest (although progress 
has been made for two dimensional systems [Henrich et 
al. (2009), Kriiger et al. (2011)]). In [Brader et al. 

(2009) ] a simplified 'schematic' constitutive model was 
proposed, which aims to capture the essential physics of 
the wavector dependent theory, while remaining numeri- 
cally tractable. Applications so far have been to steady- 
state flows, step strain and dynamic yielding [Brader et 
al. (2009)], as well as oscillatory shear [Brader et al. 

(2010) ]. 

Both the full [Brader et al. (2008)] and schematic 
[Brader et al. (2009)] mode-coupling theories predict 
an idealized glass transition at sufficiently high coupling 
strength, characterized by an infinitely slow structural 



; 1 I. INTRODUCTION 

O GoUoidal dispersions display a broad range of nontriv- 

ial rheological response to externally applied flow. Even 
the simplest systems of purely repulsive spherical col- 
loids exhibit a rate dependent viscosity in steady state 
I flows, yielding and complex time-dependent phenomena, 
such as thixotropy and ageing [Brader (2010), Mewis and 
Ch Wagner (2009)]. Understanding the emergence of these 
collective dynamical phenomena from the underlying in- 
terparticle interactions poses a challenge to nonequilib- 
rium statistical mechanics and the fundamental mecha- 
nisms involved are only beginning to be understood. The- 
oretical advances have largely been made hand-in-hand 
with improved simulation techniques [Banchio and Brady 
(2003)] and modern experimental developments, combin- 
ing confocal microscopy or magnetic resonance imaging 
with classical rheological measurements [Besseling et al. 
(2010), Frank et al. (2003)]. 

Despite considerable progress, a comprehensive con- 
stitutive theory, capable of capturing the full range of 
response, remains to be found. Existing approaches are 
tailored to capture the physics of interest within partic- 
ular ranges of the system parameters (e.g. density, tem- 
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$H perature) but fail to provide the desired global frame- 
work. Moreover, the vast majority of studies have con- 
centrated on the specific, albeit important, case of simple 
shear flow. Such scalar constitutive theories, relating the 
shear stress to the shear strain and/or strain-rate, pro- 
vide important information regarding the competition of 
timescales underlying the rheological response, but do 
not acknowledge the true three dimensional character 
of experimental flows. Tensorial constitutive equations 
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relaxation time Tq,. Ageing dynamics are neglected. An 
important prediction of the approach is that application 
of any steady strain-rate leads to fluidization of the ar- 
rested microstructure, with a structural relaxation time 
determined by the characteristic rate of flow ~ 7~^- 

In recent experiments on various soft glassy materials, 
Ovarlez et al. have indicated that when a dominant, flu- 
idizing shear flow is imposed, then the sample responds 
as a liquid to an additional perturbing shear flow, regard- 
less of the spatial direction in which this perturbation is 
applied. These flndings imply that once the yield stress 
has been overcome by the dominant shear flow, arrested 
states of soft matter become simultaneously fluidized in 
all spatial directions. In particular, the low shear vis- 
cosity in a direction orthogonal to the primary flow is 
determined by the primary flow rate. The rheometer em- 
ployed in [Ovarlez (2010)] consisted of two parallel discs 
which enabled the simultaneous application of rotational 
and squeeze shear flow, with independent control over 
the two different shear rates. Although this set-up in- 
deed provides a useful way to study superposed shear 
flows of differing rate, it does not provide a mean to test 
the three-dimensional yield surface, as claimed in [Ovar- 
lez (2010)]. A true exploration of the yield surface poses 
a considerable challenge to experiment and requires a pa- 
rameterization of the velocity gradient tensor which can 
incorporate the entire family of homogeneous flows, in- 
cluding both extension and shear as special cases. The 
superposition of two shear flows is yet another shear flow 
and does not enable the entire space of homogeneous ve- 
locity gradients to be explored. 

In the present work we will employ the constitutive 
theory of Brader et al. (2009) to investigate the response 
of a generic colloidal glass to a 'mixed' flow described 
mathematically by the linear superposition of two in- 
dependently controllable velocity gradient tensors. Nu- 
merical results will be presented for the special case in 
which simple shear is combined with uniaxial compres- 
sion. Despite the fact that we employ a combination of 
compression and shear, as opposed to the superposition 
of two shear flows, our theoretical results are broadly con- 
sistent with the experimental flndings of Ovarlez et al. 
(2010) regarding the response of shear fluidized glasses. 
In particular, our calculations reveal clearly the relevant 
timescales dictating the three dimensional response of the 
system. Following this specific application, we proceed to 
extend our description to treat more general mixed flows. 

The paper will be organized as follows: In Sec. |ll] 
we will introduce the deformation measures required to 
describe flow in three dimensions and summarize the 

(2009). 



problem of superposing a mechanical perturbation onto 
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schematic model of Brader et al. (2009). In Sec. 
will consider the application of our constitutive model to 
a specific mixed flow, namely a combination of uniaxial 
compression and simple shear. In Sec. |IV| we will present 
numerical results for the flow curves and low shear vis- 
cosity for the aforementioned flow combination. In Sec. 
[V| we will perform a perturbation analysis of our consti- 
tutive equation which enables us to address the general 



a dominant flow. Finally, in Sec. VI we will discuss the 
significance of our results and give concluding remarks. 



II. THE SCHEMATIC MODEL 

A. Continuum tensors 

Spatially homogeneous deformations are encoded in 
the spatially translationally invariant deformation ten- 
sor E{t,t'). Any given vector r{t') at time t' may be 
transformed into a new vector r{t) at later time t using 
the linear relation 



r{t) ^ E{t,t') ■ r{t') 



(1) 



where E^^p = dr^^it) / drpit') [Brader (2010)]. Calculat- 
ing the time derivative of the deformation tensor E and 
using the chain rule for derivatives yields an equation of 
motion for the deformation tensor 



dE{t,t') 
dt 



^ K{t) ■ E{t,t') , 



(2) 



where k — Vv is the velocity gradient tensor with com- 
ponents {yv)ap — dia/drp. In the present work we 
will assume incompressibility, which may be expressed 
by the condition Tr k = or, equivalently, Tr = 1 
(volume is conserved). If the deformation rate is con- 
stant in time, then the velocity gradient matrix k loses 
its time dependence {Kit) — > k) and the deformation 
tensor E becomes a function of the time difference alone 
{E{t, t') -> E{t - t')). The formal solution of Eq.([2]) for 
such steady flows is thus given by 



E{t) 



(3) 



The deformation tensor contains information about both 
the stretching and rotation of material lines (vectors em- 
bedded in the material). A more useful measure of strain 
is the Finger tensor B{t,t'), which is defined for steady 
fiows by 



B{t) = E{t) ■ E^{t). 



(4) 



The Finger tensor is invariant with respect to physically 
irrelevant solid body rotations of the material sample and 
occurs naturally in many constitutive models (e.g. the 
Doi-Edwards model of polymer melts [Doi and Edwards 
(1989)]). 



B. Schematic mode-coupling equations 

The schematic model developed in [Brader et al. 
(2009)] expresses the deviatoric stress tensor in integral 
form 



,t(<) = 



dt' 



d 



- —B{t,t')\G{t,t') 



(5) 
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An equation of the form ^ has been derived from 
first principles [Brader et al. (2008)], starting from the 
A''-particle Smoluchowski equation and applying mode- 
coupling approximations to a formally exact generalized 
Green-Kubo relation for the stress tensor. In [Brader et 
al. (2009)] the theory was simplified to ^ by assum- 
ing spatial isotropy of the modulus G{t, t'). The physical 
content of Eq.([5]) is that, in order to calculate the stress 
at the present time, increments of an appropriate, mate- 
rial objective strain measure (the Finger tensor) are inte- 
grated over the flow history, each weighted with a 'fading 
memory'. Approximating G(t,t') by an exponential re- 
covers the well-known Lodge equation [Larson (1988)], 
which is just the integral form of the upper-convected 
Maxwell model. However, Eq.([5]) differs from the simple 
Lodge equation in that, (i) the modulus G is generally not 
time translationally invariant, due time-dependent varia- 
tion of the flow in the time interval between t and t' , (ii) 
the memory does not decay exponentially to zero, but 
displays the two-step relaxation characteristic of dense 
colloidal dispersions. 

Within the wavevector dependent approach of [Brader 
et al. (2008)] the autocorrelation function of stress fluc- 
tuations is assumed to relax in the same way as the den- 
sity fluctuations. This leads to an approximation for the 
nonlinear modulus G, given by a weighted k-integral over 
a bilinear function of density correlators at two differ- 
ent (but coupled) wavevectors. The schematic model re- 
places this with the simpler form 



G{t,t') = i^„^^t,t') 



(6) 



where ^{t,t') is a single mode transient density corre- 
lator (normalized to $(t, t) = 1) and Vcr is a parameter 
measuring the strength of stress fluctuations. 

The dynamics of the single mode density correlator are 
determined by a nonlinear integro-differential equation 



$(t,to)+r|$(<,io)+ J dt'm(i,t',to)$(t',<n)| =0, 

. (7) 

where F is an initial decay rate, the inverse of which sets 
our basic unit of time. The function m(t, t' , to) is a three- 
time memory-kernel which depends upon the strain accu- 
mulated between its time arguments and describes how 
this competes with the slow structural relaxation arising 
from the colloidal interactions. The memory kernel is 
given by 



m{t,t',to) = hi{t,to)h2{t,t') 



i^i<i>{t,t') + i^2^^{t,t') 



(8) 

The dependence of the memory upon ^(t, t') is taken 
from the F12 model developed by Gotze [Gotze (2008)]. 
The coupling constants are given by vi = 2(-\/2 — 1) -|- 
e/ (\/2 — 1) and V2 ~ 2, where e is a parameter expressing 
the distance to the glass transition. The system is fluid 
for e < and in a glassy state for e > 0. 



The hi entering ([8]) are decaying functions of the accu- 
mulated strain. For simplicity we assume hi = h2 = h. 
To allow consideration of any kind of flow (not only 
shear), the function h is taken to depend upon the two 
invariants Ii and I2 of the Finger tensor 



hit,to) = 



7c 



72,+ iyhit,to) + {l-v)l2it,to)-3 



(9) 



where a mixing parameter (0 < < 1) and a cross-over 
strain parameter (7cr) have been introduced [Brader et al. 
(2009)]. The scalars h = Tr{B) and h = Tr{B^^) are 
the trace of the Finger tensor and its inverse, respectively. 
In principle, the time evolution of the density correlator 
^{t,to) and thus, via Eqs.([5| and_j6]), the stress tensor, 
can be calculated by solving Eq.(l7| numerically for any 
given velocity gradient tensor k. 

The model outlined above contains a set of flve inde- 
pendent parameters {v„,T, e, ^cr)- The least important 
of these is which determines the relative influence of 
the invariant Ii with respect to I2 in determining the 
strain induced decay of the memory function. However, 
numerical results prove to be extremely insensitive to the 
value of V, at least for all flows to which the schematic 
model has so far been applied. 

Trivial scaling of stress and time scales is provided by 
the parameters Va and F. A statistical mechanical cal- 
culation of the dynamics of N colloids (in the absence 
of hydrodynamic interactions) identifies the modulus G 
as the autocorrelation function of stress fluctuations. v„ 
therefore determines the initial value of the modulus and, 
via ([5]), sets the overall stress scale. The reciprocal of the 
initial decay rate F^^ simply acts as the fundamental 
timescale. For the purpose of our theoretical investiga- 
tions both Va and F can, without loss of generality, be 
set equal to unity. The theoretical results thus generated 
can then be fit to experimental data by scaling stress and 
time (or frequency) with alternative values for these two 
parameters [Brader et al. (2010)] 

The two most important parameters in the model are 
7cr and e. The cross-over strain 7cr sets the strain value 
at which elastic response gives way to viscous flow. For 
example, in experiments considering the shear stress re- 
sponse of dense colloidal systems to the onset of steady 
shear flow, 7cr can be identified from the peak of the 
overshoot on the stress-strain curve. The parameter e 
characterizes the thermodynamic state point of the sys- 
tem relative to the glass transition and serves as proxy for 
the true thermodynamic parameters of the physical sys- 
tem (volume fraction, temperature etc.). For example, in 
a simple system of hard-sphere colloids of volume frac- 
tion (j) one can identify e ^ (^ — (/)g)/(/)g, where 4>^ is the 
volume fraction at the glass transition. For more com- 
plicated systems e can be regarded as a general coupling 
parameter which, in the absence of flow, yields fluid-like 
behaviour for e < and amorphous solid-like response 
for e > 0. 
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FIG. 1. The flow curves of a glassy state (e > 0) for various 
values of the compressional rate 7c. The dashed line is a 
Newtonian viscous law. For % = the 7s — >■ limit of the 
flow curve identifies the dynamic yield stress [Brader et al. 
(2009)]. 



III. MIXED SHEAR AND COMPRESSIONAL 
FLOWS 

With the constitutive relation ([5]), we are in a posi- 
tion to determine the rheological behaviour of a colloidal 
glass undergoing any type of homogeneous deformation. 
In [Ovarlez et al. (2010)], Ovarlez et al. considered var- 
ious soft glassy materials loaded between two parallel 
discs. Each sample was simultaneously sheared by ro- 
tating the upper disc about its axis at a given angular 
velocity and squeezed by lowering the height of the up- 
per disc at a given rate. By independently varying the 
rotation and compression rates the stress could be deter- 
mined as a function of one of the rates, for a fixed value 
of the other. In these experiments, the rotation of the 
upper plate induces a shear flow in the (p direction (in 
cylindrical coordinates), the rate of which increases lin- 
early with radial distance from the axis of rotation. As 
a consequence of the stick boundary conditions the com- 
pression of the sample leads to an inhomogeneous shear 
flow in the f direction (somewhat akin to a Poiseuille 
flow) with a maximum shear rate at the boundaries and 
zero shear rate in the plane equidistant between the two 
plates. 

The experiments of Ovarlez et al. (2010) were per- 
formed in a curvilinear geometry using a flow protocol 
which induces an inhomogeneous velocity gradient ten- 
sor. In principle, spatial variations of the velocity gradi- 
ent could be treated within the present theoretical frame- 
work by assuming that the constitutive relations remain 
valid locally and enforcing the local stress balance appro- 
priate to the geometry of the rheometer under consider- 
ation. In addition to the increased numerical resources 
required for such an investigation, the local application 



of our constitutive equation would represent a further 
approximation, over and above those already underlying 
the schematic model. The main conceptual point emerg- 
ing from the experimental studies of Ovarlez et al. (2010) 
is that if a primary flow restores ergodicity and fluidizes 
the glass, then the response to the secondary flow is also 
fluid like. Spatial inhomogeneity of one or both flows is 
merely a complicating factor. We thus choose to focus on 
a more idealized homogeneous flow which is convenient 
for numerical implementation, but nevertheless captures 
the salient features of the experiment in a minimal way. 

The homogeneous flow we choose to implement is a 
superposition of simple shear and uniaxial compressional 
flow. We anticipate that the key physical mechanism 
at work in fluidized systems under superposed flow is 
the competition between the two imposed relaxation 
timescales. As the superposition of two shear flows is 
itself another shear flow, the experiments of Ovarlez et 
al. (2010) leave open the possibility that the observed 
phenomena could be a special feature of shear. For this 
reason we chose to implement the mathematically more 
general case of superposed extension and shear, for which 
the geometrical coupling of the flows is more involved. 

Working in a cartesian coordinate system our flow is 
specifled by 



(10) 
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The shear and compressional flows are represented by the 
following matrices 



(11) 



where % and 7c are the shear and compression rates, re- 
spectively. Our choice of flow thus differs from those of 
Ovarlez in two respects, (i) both 7^ and 7c are trans- 
lationally invariant and, (ii) we superpose shear with 
genuine elongation, as opposed to superposing two shear 
flows. We consider the flow ( [To| as a thought experiment 
intended to highlight the fundamental physical mecha- 
nism of fluidization in a simple and transparent fashion. 



A direct experimental realization of (10 1 is not feasible, 
as this would require a rheometer with stick boundary 
conditions for generating the shear flow, but slip bound- 
aries for the compressional flow. As we will see below, 
our assumptions do not seem to lead to qualitative dif- 
ferences between our theoretical findings and the exper- 
imental results and simplify considerably the theoretical 
calculations. 

Eq.([3]) enables calculation of the deformation tensor 
E{t) for our mixed flow. The non-zero elements are given 
by 

p7ct/2 



E^, 



E, 



-7ct 



2% 
37c 



(12) 
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FIG. 2. The low shear viscosity r/ scaled by the yield stress 
(jyicid as a function of compression rate 7c (glassy states with 
£ = 10~^,10~^ and 10~^). The continuous line is a power- 
law fit to the numerical data points for e = 10~^ over the 
range 7c = 10~^^ — 10""* and yields an exponent of -1. The 
e-dependent deviations apparent for 7c > 10~* indicate that 
short-time relaxation processes are becoming relevant. 



Employing Eq.Q yields the Finger tensor 
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with inverse given by 
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(13) 



(14) 



The invariants required for the memory function prefac- 
tors ([9| are thus 

hit) ^ 2e^-' + e-^^^' + Ely , (15) 

hit) = 2e-^=* + 6^^=* + Elye^^' . (16) 

Finally, we need to calculate the time derivative of the 
Finger tensor B{t). In Sec. |IV| we will present results for 
the shear stress (7 xy SI'S 01 function of treating 7c as a 
parameter. Inspection of ([s]) shows that we require only 
the xy component of the Finger tensor time derivative 



dB,y{t) 
dt 



Iip-27, 

3 



(17) 



Substituting (17) into ([s]) and assuming time trans- 
lational invariance (as appropriate for the steady flows 
under consideration) we obtain our final expression 



'xy 



= / dt 



Iip-27, 

3 



:^^<P\t). (18) 



The xy component of the shear stress tensor is now com- 
pletely cha ract erized. When numerically evaluating the 
integral in ( 18 ) we find that truncation at t ^ (is+ic)^^ 



provides accurate results. We note that, in an analogous 
way, all other components of the shear stress c(i) can 
be calculated, which is useful if one is interested for ex- 
ample in the first and second normal stress differences, 
iVi = a^x - cTyy and = cfyy - (^zz , respectively. 



IV. NUMERICAL RESULTS 

In Figjljwe show fiow curves generated from numerical 
solution of Eqs.( 7][9 ) and ( 18 ). For each curve we set the 
compressional rate equal to a fixed value, in effect treat- 
ing 7c as a parameter, and plot the shear stress axy as a 
function of The model parameters used to generate 
these data are as follows: (F — 1,1/^ — l,7cr — — 
0.5, e = 10~^). For 7c = we recover the simple shear 
flow curve which, for the glassy state under consideration, 
tends to a dynamic yield stress a^y — >■ fyieid = 0.03633 in 
the limit of vanishing shear rate. Within the theory, the 
existence of a dynamic yield stress is a direct consequence 
of the scaling of the structural relaxation time with shear 
rate, ^ 77^- The flow curves calculated at finite 7c 
differ qualitatively from that at 7c = 0. In particular, 
CTyioid is a discontinuous function of the parameter 7c, 
such that f7:cy(7s -> 0,7c = 0) 7^ (yxyiis -> 0,7c -> 0+). 
For finite 7c values the fiow curves present a Newtonian 
regime for rates 7s < 7c, followed by a shear thinning 
regime for 7^ > 7c. The existence of two regimes is quite 
intuitive: For 7s < 7c compression is the dominant, i.e. 
fastest, fiow and sets the timescale of structural relax- 
ation, whereas for 7s > 7c the shear flow dominates and 
the flow curve converges to the 7c = result. 

The above findings are in good qualitative agreement 
with the experimental results obtained in [Ovarlez et al. 
(2010)] (cf. Fig. 3 therein). In order to characterize more 
precisely the flow curves at flnite 7c we show in Figj2] 
the low shear viscosity — (T^y/'^sj scaled by the yield 
stress CTyicid, as a function of 7c for three different positive 
values of e. For 7c < 10~^ we find very good data collapse 
onto a master curve. For 7c > 10^^ clear deviations 
from universality set in, signifying that the compression 
induced structural relaxation processes are occurring on 
a timescale within the microscopic regime, for which 77 
becomes an e independent quantity (around % = 10"'* 
for the parameter set used in Figjl]). Provided that 7c < 
we find that the numerical data are well represented 



by the power-law scaling 



(19) 



with a = —1, in agreement with the experimental find- 
ings of Ovarlez et al. (2010). The constant of propor- 
tionality A is independent of e (both 77 and tXyicid vary 
in the same way with this parameter). Given the lack 
of detailed material specificity in the schematic model. 
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which point ^cff = e > 0. The inset of FigjS] shows the 
same data on a logarithmic scale. In this representation 
it becomes apparent that the data follow a power law 



0.015 



FIG. 3. Mapping a compression melted glassy state onto an 
effective fluid state obtained by matching the low shear vis- 
cosity. The effective distance to the glass transition Scft is 
shown as a function of the compressional rate 7c, for different 
values of e. The points are numerical data and the lines pro- 
vide a guide for the eye. The inset shows the same data on a 
logarithmic scale. 



we are led to believe that a = —1 is a universal ex- 
ponent, independent of both the details of the material 
under consideration and of the precise nature of the pri- 
mary and perturbing flows. Our findings suggest that any 
constitutive theory capable of describing a three dimen- 
sional dynamic yield stress ('yield stress surface' [Brader 
et al. (2009)]) will inevitably recover the scaling (19) 
with a 



-1, when applied to tackle mixed flows. In 
particular, we anticipate that the full wavevector depen- 
dent mode-coupling constitutive equation [Brader et al. 
(2008)] would predict the same scaling behaviour, al- 
though this claim remains to be confirmed by explicit cal- 
culations. Within mode-coupling-based approaches the 
value of the scaling exponent a is a natural consequence 
of the way in which strain enters the memory function 

The fiow curves presented in Fig[l]for various values of 
7c are very reminiscent of the (more familiar) fiow curves 
either measured or calculated under simple shear with 
7c = and £ < 0, i.e. states which would remain fiuid in 
the absence of fiow (see, e.g. Fuchs and Gates (2003)). 
This similarity suggests that it may be possible to map, 
at least approximately, the shear response of a steadily 
compressed, glassy system with 7c 7^ and e > onto 
an uncompressed, fluid system, t-c = 0, at some effective, 
negative value of the separation parameter Scff- One pos- 
sible way to realize such a mapping is to adjust Eeff for a 
given 7c to obtain equal values for the low shear viscosity 
of the compressed glass and effective fluid systems. The 
results of performing this procedure for three values of e 
are shown in FigjSj It should be noted that the mapping 
between Soff and 7c becomes discontinuous at 7c = at 



(20) 



Fits to our numerical data yield values for the exponent 
0.41 < /3 < 0.43. 

Within the quiescent F12 schematic model [Gotze 
(2008), Gotze (1984)], to which the present theory re- 
duces in the absence of flow, it is known that the zero 
shear viscosity exhibits a power law divergence as e ap- 
proaches the glass transition from below 



(21) 



where 5 is the same exponent as that describing the di- 
vergence of Ta at the glass transition. Note that the sym- 
bol 5 is employed here for this exponent, rather than the 
standard choice 7, in order to avoid confusion with the 
strain. When employing the Percus-Yevick approxima- 
tion to the static structure factor as input, the wavevector 
dependent mode-coupling theory predicts that for hard- 
spheres the viscosity exponent takes the value 5 — 2.46 
(identifying e as the volume fraction, relative to the tran- 
sition point) [Gotze and Sjogren (1992)]. Within the 
present schematic model we obtain 5 = 2.3 (see also foot- 
note [H]). Given this information about the divergence 
of rj in the quiescent system, the power law relation for 
the mapping ( 20 ) is already implicit in the data shown 
Figj2] Using the relations (19 1 and (21) the relation 



( |20[ ) can be deduced, where the exponent fi is given by 
/3 = a/{—5) = 0.43, which is consistent with the results 
of our numerical fits. 



V. ANALYTIC PERTURBATIVE RESULTS 

We have so far focused on the special case of mixed 
shear and compressional fiows. For any given value of 
7c < F we have shown that there exists a Newtonian 
regime in the stress response to the shear flow, provided 
is < 7c (see Figjl]) [32]. In this section, we now con- 
sider more general situations for which a second slow 
flow K2 is added to a dominant flow Ki (while keep- 
ing the requirements of incompressibility and homogene- 
ity). In the present context a sufficient condition for 
the second flow to be considered 'slow' is that 72 ^71, 
where t he char acteristic shear rates are now identified as 
7i = ^/K~^K[ fori = 1,2 (where A:B = J^ij^ij^ji)- 
In Subsections |V A]V C[ we provide perturbative consti- 
tutive equations for three different cases. In the first 
of these cases, we consider Hi and K2 as steady flows 
(without any other restriction), derivate the correspond- 
ing pertubative constitutive equation and flnally apply 
this latter to our coupled compressional and shear flows, 
in order to theoretically account for the Newtonian vis- 
cous response to K2 discussed in Sec. |IV[ and to finally 
make the connection with the phenomenological consti- 
tutive equation obtained by Ovarlez et al.. In the second 
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case, we still consider steady flows, but this time with 
the additional requirement of 'commutating' flows, i.e. 
[/ti, K2] = ti • /«2 ^ '*2 ■ '^i = 0, whereas in the third case 
K2 is time-dependent and the requirement of commutat- 
ing flows is maintained. We also illustrate these last two 
cases with instructive examples. 



A. Newtonian viscous response 

1. Perturbation expansion 

The first point to note is that the isotropic modulus 
G(t) decays on the timescale 7]"^; the slower secondary 
flow has no influence on the structural relaxation. We 
henceforth make this fact explicit in the notation for the 
modulus by writing G'(t; 71, 72) = G{t;^i) = Gi{t). For 
steady flows the constitutive equation ([5| may thus be 
simplified to 



dt 



d 



B{t) G^{t). 



(22) 



Although the modulus is essentially independent of 72, 
the Finger tensor depends nonlinear ly upon both Hi and 
K2. In order to address the case 72 <C 71 we expand the 
Finger tensor Q to first order in k,2 ■ For the mixed flow 
under consideration B{t) is given by 



B{t) 



(23) 



The desired partial linearization of ( [23| ) with respect to 
K2 is complicated by the fact that the two velocity gra- 
dient tensors do not necessarily commute. 

In order to proceed we consider the following Taylor 
expansion 



a 



da 



+ 0(«'), (24) 



where x and y are arbitrary operators independent of 
the scalar coupling parameter a. The derivative may be 
obtained using the Feynman identity [Feynman (1951)] 



^x+ay 



da 



Q = "'0 



(25) 



Applying (24) and (251 to (23 1 we obtain the leading 
order result 



(26) 



B{t) ^ Bi{t) + J ds [U{t, s) + {t, s 
where we define the following tensors 

U{t,s) = Ei{t)-Ei{-s)-K2-Ei{s)-El{t), 



(27) 
(28) 



where Ei(t) = exp(/«it). Eq.(26) is linear in K2 but 



retains all orders of the dominant flow Ki. Substitution 



of ( 26 1 into ( 22 ) thus yields a stress tensor consisting of 
d 



two terms 
(T = dt 







dt 



dt 
d_ 
dt 



Bi{t) Gi{t) 



ds 



(U{t,s) + U^{t,s)) 



Gi{t) 



CTl + 5CT 



(29) 



where cri is the stress arising purely from the dominant 
flow and 5cr is the additional contribution from the slow 
perturbation. 



The perturbative term 5a in (29) is a tensor whose 
elements depend upon time as [•] 



where Uij > 1 



is an integer. Within the schematic model the relax- 
ation time determining the decay of Gi{t) is given by 
Ta = 7cr/7ii where 7cr is the cross-over strain parameter 
entering (|9|. This decay serves to cut off the integral 
in ( |29[ ) at the upper limit t ~ Tq., with the consequence 
that the numerically largest elements of arising from 
terms with n.y = 1 are C'(?7i7cr). In a (repulsive) col- 
loidal glass any given colloid is trapped within a cage 
of nearest neighbours. The cross-over strain parameter 
7cr is related to the strain at which the cages begin to 
be broken by the external flow. Typical values for this 
dimensionless parameter from simulation or experiment 
are 7cr « 0.1 [Zausch et al. (2008)]. 



We now apply the perturbative formula ( 29 ) to the 



coupled compressional and shear flows expressed by the 
matrices (11 ), with Ki — and k,2 — Kg. Since no com- 



plete analytic expression is known for the density corre- 
lator <&, we approximate the modulus Gi{t) by an expo- 
nentially decaying function Gi{t) w Goo exp(— 7it/7cr), 
where Goo is a constant. Under this approximation, Eq. 



(29) becomes 







l+27c 





2vi 



(2 + 57c. + 2jI) 



Do 



(30) 



where r/i = dt Gi{t) is the rate dependent shear vis- 
cosity of the primary flow alone and D2 = (/«2 + 1^2)/'^ 
is the symmetric part of the velocity gradient matrix K2- 
The second term in ( 30 ) is nothing but the expression of 



a Newtonian-type viscous response to the secondary flow 
K2 (shear flow) with a viscosity mainly determined by 
the strain rate of the dominant flow 71 (compressional 
rate) through rji. This is in agreement with what we 
numerically showed in Sec. |IV| 



2. Empirical constitutive equation 

In [Ovarlez et al. (2010)], Ovarlez et al. proposed an 
empirical constitutive equation to account for the viscous 
stress measured in a number of fluidized glassy systems. 
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In the notation of the present work the proposed consti 
tutive relation is 

kd 



cr = 2 



fyicld 



d 



(31) 



where k and n are scalar parameters and d = \/2D : D 
is an invariant of the symmetric part of the total ve- 
locity gradient D = {k + k'^)/2. The isotropic vis- 
cosity appearing in square parentheses in (31) is ob- 
tained from a straightforward generalization of the fa- 
miliar scalar Hershel-Bulkley law for the shear stress, 

If we neglect the cross-over strain parameter j^r (whose 
value is already small, jcr ~ 0.1), then the second term 



in our perturbative constitutive equation ( 30 1 becomes 
Scr = 2riiD2, which is entirely consistent with the im- 
plicit second term in the empirical relation (31 ). Indeed, 



for 7i ^ 72 the generalized Hershel-Bulkley effective vis- 
cosity (o-yicid + kd^^)/d is dominated by the fastest flow 
and is effectively independent of 72 . We can thus make 
the following correspondence between viscosities appear- 



ing in the schematic (30) and empirical (31 1 constitutive 
equations 



dtGi{t) 



fyiold 



kd" 



(32) 



The linear dependence of D on the velocity gradient ten- 
sor thus enables (31 1 to be rewritten as 



(33) 



thus making explicit the connection between (30) and 

dsTl). 



B. Anisotropic viscosity 

We still consider steady flows Ki (dominant flow) and 
K2 (secondary flow), but now with the restriction of 
commutating flows, i.e. [k;i,/«2] = 0. Such flows have 
the property that the total deformation tensor can be 
formed from the product of the individual deformations, 
E{t) ~ Ei{t)E2{t). As we will see, this restriction allows 
for more tractable pertubative constitutive equations. 

With [/«i,/«2] = 0, the expression (26) then reduces to 



B{t) = Bi{t) + 2tEi{t)-D2-El{t) 



(34) 



Substitution of (34) into (22) yields the following form 
for the stress tensor 



cr = o"! + 21J1D2 



(35) 



2 / dtGi{t) 
/o 



d_ 
dt 



{tE^{t)-D2-El{t)^-D2 



The anisotropic third term in ( 35 ) is the result of a non- 
linear operator acting on the perturbing velocity gradient 
D2 and incorporates information about the symmetry 
imposed on the system by the dominant fluidizing flow. 



For incompressible isotropic fluids in the Newtonian 
regime the viscosity in any given flow can be determined 
from the shear viscosity via Trouton's rules (e.g. r^oi = 
3?7sh, where ?7oi is the elongational viscosity in uniaxial 
extension). Trouton's rules no longer hold in the present 
case, due to the presence of the third term in (351. 



In order to explicitly demonstrate the relative mag- 
nitude of the anisotropy, we consider the special case 
of perpendicular shear flow = jiSixSjy, {i^2)ij = 

i2SizSjy and again approximate the modulus by an expo- 
nentially decaying function Gi{t) w Goo exp(— 7it/7cr). 



Under this simplifying assumption Eq. ( 35 ) becomes 



(T = (Tl+ 277172 




1/2 
1/2 








7cr \ 


1+ 277172 1 





s) 










(36) 



The small off-diagonal elements which appear in the third 
term of ( 36 ) are generated by the coupling between pri- 



mary and perturbing flows and may be viewed as a 
correction, at around the 10% level, to the dominant 
isotropic viscosity 771 . The appearance of these additional 
contributions to the viscous stress can be attributed to 
the normal stress differences generated by the primary 
flow. Constitutive theories with vanishing normal stress 
differences will always predict an isotropic viscous re- 
sponse to perturbing flows. 



It is interesting to note that the stress tensor ( 35 ) may 



be formally expressed in terms of an anisotropic viscosity 

{(T)ij = (cri)jj- + 2{T])ijki{D2)ki, (37) 

where the fourth rank tensor r) with components {ri)ijki is 
given in terms of the shear modulus and the deformation 
gradient of the dominant flow 



(r7),,fei = 2/ dtGiit) 
Jo 



t{E,{t)U{E,{t))^ 



3 1 



(38) 



If the dominant flow is switched off, then Ei (t) ~ 1 and 
(38) reduces to the familiar isotropic viscosity rjijki = 
'^Vo^ikSji, where 770 is the zero shear viscosity (infinite 
for glassy states with e > 0). Given that the dominant 
flow fixes the anisotropy of the system it is not surprising 
that the viscosity experienced by the perturbing flow is 
a tensorial quantity. 

Finally, we note that the presence of anisotropy pre- 
vents a general three dimensional mapping of a flow flu- 
idized glass onto an effective fluid state with e < 0, as 
performed for the special case of mixed compressional 
and shear flow in Sec. IIVI 



C. Superposition spectroscopy 

A special case of mixed flow which has received some 
attention in the rheological literature is small amplitude 
oscillation superposed onto steady shear. Largely due 
to constraints imposed by the available apparatus, the 
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FIG. 4. The othogonal superposition storage (G'l, filled sym- 
bols) and loss (G", open symbols) moduli as a function of 
frequency for a fiuid state. The (black) circles show the stan- 
dard linear response moduli calculated at a vanishing steady 
shear rate 71 = 0. The (red) squares and (blue) triangles show 
the moduli calculated at finite steady shear rates, 71 = 10~^ 
and 10"'', respectively. The arrows indicate the character- 
istic frequency 27r7i below which the structural relaxation 
is dominated by the steady shear (which coincides approxi- 
mately with the crossing point of the moduli) . Parameter set 
(r = 1,1/,, = 100, 7„ = 1,!/= 1,£ = -5 X 10^3). 



majority of the experimental works have involved paral- 
lel shearing flows using either cone-plate [Booij (1966), 
Osaka et al. (1965)] or Couette [Vermant et al. (1998)] 
rheometers. The measured viscoelastic parallel superpo- 
sition moduli depend upon both the microstructure un- 
der steady shear and its evolution with changes in shear 
rate. 

A more informative, albeit harder to realize, mixed 
flow consists of oscillatory shear superposed perpendicu- 
lar to the main flow direction. The orthogonality of the 
flows makes possible a mechanical spectroscopy of flow- 
ing systems which can probe flow induced changes in the 
microstructure and which may be used as a useful test 
of constitutive equations [De Cleyn and Mewis (1987), 
Kwon and Leonov (1993), Leonov et al. (1976), Tanner 
and Simmons (1967)]. For details on the experimental 
realization of such orthogonal flows we refer the reader 
to Vermant et al. (1997). 

As a specific example of orthogonal oscillation we con- 
sider the mixed flow 



(39) 
(40) 



where uj is the angular frequency and 70 is the amplitude 
of the oscillatory strain (assumed to be small) . The time- 
dependence of the perturbing flow K2 requires us to use 
a time-ordered exponential to express the corresponding 



deformation tensor, namely 

E,{t,t')^J_^''''^^'^'^ , 
where the exponential is defined according to 



(41) 



n,dsA{s) 



^1+ dsA{s)+ dsi / ds2 A{si)A{s2)+. 
Jt' Jt' Jt' 



(42) 

We also require the time-dependent expression for the 
stress tensor given by ([s]). After linearization with re- 
spect to K2 and making use of [ki, K2] = 0, we obtain a 
formula for the stress tensor analogous to ( 35 ) , 



cr(t)=cri+ f dt'Gi{t,t'-^i) 

■/ — 00 

A \^E,{t - t')-(^fds 2 r>2(.s)) -E^^it t')] } 



(43) 



Substituting ( 39 ) and ( 40 1 into ( 43 1 and making use of 



standard trigonometric addition formulas yields 



' zy 



(t) =7oGi(w; 71) sin(wt) + 7oG"(w; 71) cos(wi). 



(44) 



where the orthogonal superposition moduh are given by 
G[{uj;ji)=uj dt'sin(wt')Gi(t';7i), (45) 



G'/(w;7i)=w/ dt'cos{ujt')Gi{t';ji). (46) 
Jo 



In Eqs.( 44p6 1 we have made explicit the dependence of 
the moduli upon the steady shear rate 71. The applica- 
tion of oscillations perpendicular to the fiow thus enable 
the modulus under steady shear to be investigated and 
provide information about the shear induced relaxation 
of stress fluctuations. We note that identical moduli (45) 
and (46) would be obtained had we chosen the alter- 
native perturbing flow (/t2)ij = Jo'-^ cos{ujt)SixSjz and 
determined the stress component axz- 

In Fig|4]we show the orthogonal superposition moduli 
as a function of frequency for three different values of the 
steady shear rate 71 . For 71 = we recover the standard 
linear response moduli, for which the viscous loss domi- 
nates the elastic storage for frequencies less than 27r/rQ. 
For the fluid state considered Tq is finite (e < 0). As 
the steady shear rate is increased, relaxation processes 
with rates less than 71 are suppressed and the point at 
which the storage and loss moduli cross moves to higher 
frequency. These findings are consistent with the exper- 
imental results of both Booij (1966) and Vermant et al. 
(1998). The underlying physics here is essentially the 
same as that leading to the shift of the Newtonian regime 
shown in Figjl] At low frequencies the orthogonal super- 
position moduli retain the same frequency scaling as in 
the familiar unsheared situation, namely G'l '-^ oj^ and 
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G'( ~ w. We note also that the Kramers-Kronig rela- 
tions remain valid for finite values of 71 . 

In addition to the speeding up of structural relaxation 
induced by the steady shear, the loss modulus also dis- 
plays a more pronounced a-peak compared to the un- 
sheared function. This feature is related to the functional 
form of the a-decay of the transient density correlator. 
In the absence of shear ^{t) decays as a stretched expo- 
nential, whereas under shear the final decay is closer to 
pure exponential. However, it is likely that more accu- 
rate (i.e. beyond schematic) orthogonal moduli, obtained 
either from experiment or more detailed microscopic cal- 
culations/simulations, would differ qualitatively in the 
region of the a-peak. There is accumulating evidence 
[Zausch et al. (2008)] that Gi(t;7i) becomes negative 
at long times and, as this feature is not captured by the 
simple schematic model employed here, differences in the 
Fourier transformed quantity around w = 27r7i may be 
anticipated. The negative tail of Gi(f;7i) is related to 
the existence of a maximum in the shear stress as a func- 
tion of time following the onset of steady shear [Zausch 
et al. (2008)]. The 'stress overshoot' is present in the full 
wavevector dependent mode-coupling equations [Brader 
et al. (2008)] but gets lost in simplifying the theory to 
the schematic level. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have demonstrated that the MCT- 
based schematic model of Brader et al. (2009) can quali- 
tatively account for the experimental results on three di- 
mensional fiow of soft glassy materials reported in [Ovar- 
lez et al. (2010)]. In particular, the competition of 
timescales which arises from applying fiows of differing 
rate appears to be correctly incorporated into the model. 
The main outcome of our analysis is that the viscous re- 
sponse to a perturbing secondary flow is dominated by 
the primary flow rate. Although subtle anisotropic cor- 
rections to this picture do emerge from our equations, it 
remains to be seen whether these have significant con- 
sequences for experiments in any particular rheometer 
geometry. 

A key feature of the mode-coupling constitutive theory 
is that it captures the transition from an ergodic fiuid to 
an arrested glass as a function of the coupling strength. 
The present study demonstrates that the theory quali- 
tatively accounts for experimental data on mechanically 
fiuidized glassy systems in three dimensional situations. 
What remains to be established is whether the experi- 
mental yield surface of a colloidal glass agrees with the 
(almost) von Mises form [Hill (1971)] predicted by the 
schematic model [Brader et al. (2009)]. A true measure- 
ment of the yield surface would require a rheometer which 
enables parameterization of the entire family of velocity 
gradients, incorporating both uniaxial and planar exten- 
sional fiows. This has not yet been achieved. Given the 
very different mechanisms underlying plastic flow in col- 



loidal glasses and metals (for which the von Mises yield 
surface was originally proposed) direct measurements of 
the yield surface could prove very informative. 

The good qualitative agreement of our theory with the 
experimental results of Ovarlez et al. (2010) on fiuidized 
glasses is perhaps all the more surprising when recalling 
that the theory is constructed specifically for dispersions 
of spherical colloidal particles (without hydrodynamic in- 
teractions), whereas the experiments were performed on 
large aspect ratio Bentonite clay, a Carbopol gel and an 
emulsion. The consistent phenomenology presented by 
these disparate systems would seem to indicate that the 
sufficient elements required for a successful theory are (i) 
a well-founded geometrical structure (in the sense that its 
tensorial structure is appropriate) , (ii) correct incorpora- 
tion of fiow induced relaxation rates. The specific nature 
of the interparticle interactions does not seem to be of 
particular importance, although we note that certain in- 
teraction potentials may be more susceptible to inhomo- 
geneous flow (e.g. shear banding instabilities) than oth- 
ers. The presence of a spatially varying velocity gradient 
tensor k would conflict with the assumption of transla- 
tional invariance underlying our constitutive equation. 

For the case of superposed compression and shear flow 
we have found that the viscosity felt by the perturbing 
shear flow is given by 277i/((2 + 5jcr + '^Icr))^ which for 
the typical value jcr = 0.1 is around 20% less than the pri- 
mary viscosity rji. In [Ovarlez et al. (2010)] the sedimen- 
tation velocity of a sphere falling in the vorticity direction 
of a shear fluidized glass was observed to be a factor of 
1.4 larger than one would expect from a sphere falling 
through a fluid of viscosity rji. Ovarlez et al. attributed 
this to hydrodynamic interactions between sedimenting 
particles. Although the flow around a falling sphere in 
shear flow is more complicated than the fiows consid- 
ered in the present work, it is nevertheless tempting to 
speculate that enhanced sedimentation velocity could be 
connected to a reduced effective viscosity, as occurs in 
( |30| , arising from a nontrivial coupling of the superposed 
flows. We leave a detailed application of our constitutive 
equation to the problem of sedimentation under shear to 
future work. 

An aspect of the present work which may warrant fur- 
ther investigation is the possible analogy between sys- 
tems with isotropic interparticle interactions, upon which 
anisotropy is imposed by external mechanical force fields, 
and intrinsically anisotropic materials such as liquid crys- 
tals [de Gennes and Prost (1993), Larson (1988)]. The 
theory of anisotropic fiuids has a long history, begin- 
ning essentially with the work of Oseen in the 1930s [Os- 
een (1933)] and developed through the work of Ericksen 
(1959) and Leslie (1968). In all of these theoretical devel- 
opments the anisotropy of the viscous response originates 
from the underlying anisotropy of the constituents; usu- 
ally oriented polymers or rod-like particles with liquid 
crystalline order. Take the nematic phase as an example. 
Within a continuum mechanics description the orienta- 
tional order is characterized by the director n. In certain 
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systems the director may be held fixed by the apphcation 
of a suitably strong external field, in others it interacts 
with the imposed flow in a more complicated way. Either 
way, the presence of a preferred direction in the sample 
gives rise to an anisotropic viscous response (character- 
ized e.g. by the five scalar Leslie viscosity coefficients). 
In the present case the anisotropy of the perturbation 
response is determined by the geometry of the primary 
flow. It may thus be anticipated that the eigenvectors of 
the primary deformation may play a role in the present 
theory analogous to that of n in the dynamics of nemat- 
ics. 

Throughout the present work we have focussed on the 
response of a glass which has already been fiuidized by a 
primary flow of constant rate. However, within the same 
formalism we can also consider the predictions of our con- 
stitutive equation for the elastic response of a colloidal 
glass which has been pre-strained by a primary defor- 
mation at some point in the past (see the Appendix for 
more details on this point). For example, a colloidal glass 
subject to a shear strain below the yield strain may rea- 
sonably be expected to possess an anisotropic elastic re- 



sponse to an additional small perturbing strain (see ( 50 ) 



for verification of this assertion) . Taking this idea a step 
further, it would be of interest to investigate the nature 
of the yield stress surface in such pre-strained glasses. 
Although both the topology of the surface and its in- 
variance with respect to hydrostatic pressure will remain 
unchanged by pre-straining, significant deviations from 
circularity (over and above those already arising from 
normal stresses) can be envisaged. In reality these devi- 
ations may well be nonstationary, decaying away as the 
sample ages, but such subtle dynamic effects are beyond 
current formulations of the mode-coupling theory. 
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Anisotropic elasticity 

Integrating ([s]) by parts yields the following form for 
the schematic constitutive equation 



a(t) ^B(t, -oo)G(i, -oo) - ly^l 



+ y* dt' B{t,t')-^,G{t,t'). 



(47) 



For glassy states the modulus relaxes to a plateau value 
for long times and the contribution of the integral term 
to the overall numerical value of the stress is limited to a 
negligible integration over the j3 relaxation (beyond this 
time the time derivative vanishes). The integral term 
may therefore be neglected to a good level of approxima- 
tion. 

We first consider the situation when the system is 
subjected to a mixed strain field e = ei + e2, where 
Ci = (ViTj -l-r^Vi) is the infinitessimal strain due to flow 
i. Both of these strains are sufficiently small that the 
system remains in the purely elastic regime. The stress 
for long times after the application of the two strains is 
given by a simple linear superposition of the two elastic 
responses 



<T — 2Gooei -|- 2Gas^2, 



(48) 



as is expected from the linear theory of isotropic elasticity 
[Landau and Lifshitz (1986)]. Note that in (48) we have 
suppressed an irrelevant isotropic contribution —v^l (the 
system is incompressible). 

We next consider the situation whereby e2 remains 
small but the strain due to the primary flow is allowed 
to be sufficiently large that some plastic rearrangements 
are induced. We nevertheless require that the total strain 
must still remain sufficiently small that the yield stress is 
not exceeded and the system remains solid. The analy- 
sis of this case neccessitates use of the partial expansion 
(34). We assume that the primary strain has been ap- 



plied at some time in the distant past and that all plastic 
rearrangements have ceased by the time we apply the per- 
turbing strain 62- The Finger tensor Bi{t) = Bi is thus 
independent of time at the present time t. For times t 
following application of £2 



T — (T 1 -f 2Goo Ei-e2- El . 



(49) 



By analogy with the situation considered in Sec. VB 



Eq. ( 49 ) can be expressed in terms of a fourth rank stiff- 



ness tensor 



{cr),j = {a- 



(50) 



where {C)ijki = 2Goo{Ei)ik{Ei).ji. We note that the 
plateau value of the modulus Goo may differ from its qui- 
escent value as a result of plastic deformation occurring 
during or after application of the primary strain. 
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